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This chapter discusses the dynamical properties of the Kuiper belt population. Then, it 
focuses on the characteristics of the Kuiper belt that cannot be explained by its evolution 
in the framework of the current solar system. We review models of primordial solar system 
evolution that have been proposed to reproduce the Kuiper belt features, outlining advantages 
and problems of each of them. 



1. Introduction 

Since its discovery in 1992, the Kuiper belt has slowly 
revealed a stunningly complex dynamical structure. This 
structure has been a gold mine for those of us interested in 
planet formation because it provides vital clues about this 
process. This chapter is a review of the current state of 
knowledge about these issues. It is divided in two parts. 
The first part (section|2]l is devoted to the description of the 
current dynamics in the Kuiper belt. This will be used in 
section |3]to highlight the properties of the Kuiper belt pop- 
ulation that cannot be explained by the current dynamical 
processes, but need to be understood in the framework of a 
scenario of primordial evolution of the outer Solar System. 
In the second part, we will review the models that have been 
proposed so far to explain the various puzzling properties 
of the Kuiper belt. More precisely, section |4] will focus on 
the origin of the outer edge of the belt; section |5] will de- 
scribe the effects of the migration of Neptune on the orbital 
structure of the Kuiper belt objects (KBOs) and section |6] 
will discuss the origin of the mass deficit of trans-Neptunian 
population. In section [7] we will present the consequences 
on the Kuiper belt of a model of the outer Solar System 
evolution that has been recently proposed to explain the or- 
bital architecture of the planets, the Trojan populations of 
both Jupiter and Neptune, and the origin of the Late Heavy 
Bombardment of the terrestrial planets. A general discus- 
sion on the current state-of-the-art in Kuiper belt modelling 
will close the paper in section [8] 

2. Current dynamics in tlie Kuiper belt 

Fig. [T] shows a map of the dynamical lifetime of trans- 
Neptunian bodies as a function of their initial semi-major 
axis and eccentricity, for an inclination of 1° and with their 



orbital ellipses oriented randomly (Duncan et ai, 1995). 
Additional maps, referring to different choices of the ini- 
tial inclination or different projections on orbital element 
space can be found in Duncan et ai, 1995 and Kuchner et 
ai, 2002. These maps have been computed numerically, 
by simulating the evolution of massless particles from their 
initial conditions, under the gravitational perturbations of 
the giant planets. The planets were assumed to be on their 
current orbits throughout the integrations. Each particle 
was followed until it suffered a close encounter with Nep- 
tune. Objects encountering Neptune would then evolve in 
the Scattered disk (see chapter by Gomes et ai). 

In fig. [T] the colored strips indicate the timespan required 
for a particle to encounter Neptune, as a function of its ini- 
tial semi-major axis and eccentricity. Strips that are colored 
yellow represent objects that survive for the length of the 
simulation, 4 x 10^ years (the approximate age of the Solar 
System) without encountering the planet. The figure also 
reports the orbital elements of the known Kuiper belt ob- 
jects. Green dots refer to bodies with inclination i < 4°, 
consistent with the low inclination at which the stability 
map has been computed. Magenta dots refer to objects with 
larger inclination and are plotted only for completeness. 

As the figure shows, the Kuiper belt has a complex dy- 
namical structure, although some general trends can be 
easily identified. If we denoted the perihelion distance of 
an orbit by q, and we note that q — a{l — e), where a is the 
semi-major axis and e is the eccentricity, Fig.[T]shows that 
most objects with q < 35 AU (in the region a < 40 AU) or 
q < 37-38 AU (in the region with 42 < a < 50 AU) are 
unstable. This is due to the fact that they pass sufficiently 
close to Neptune to be destabilized. It may be surprising 
that Neptune can destabilize objects passing at a distance of 
5-8 AU, which corresponds to ~ 10 times the radius of Nep- 
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Fig. 1. — The dynamical lifetime for massless particles 
in the Kuiper belt derived from 4 billion year integrations 
(Duncan et ai, 1995; but extended to a = 60 AU for this 
review). Each particle is represented by a narrow vertical 
strip of color, the center of which is located at the particle's 
initial eccentricity and semi-major axis (the initial orbital 
inclination for all objects was 1 degree). The color of each 
strip represents the dynamical lifetime of the particle, as re- 
ported on the scale on the right hand side. For reference, 
the locations of the important Neptune mean-motion reso- 
nances are shown in blue and two curves of constant peri- 
helion distance, q, are shown in red. The (a, e) elements of 
the Kuiper belt objects with orbits determined over 3 oppo- 
sitions are also shown. Green dots are for i < 4°, magenta 
dots otherwise. 

tune's gravitational sphere of influence or Hill radius. The 
instability, in fact, is not due to close encounters with the 
planet, but to the overlapping of its outer mean motion res- 
onances. It is well known that mean motion resonances be- 
come wider at larger eccentricity (see for instance Demon 
and Murray, 1983 and also Morbidelli, 2002) and that res- 
onance overlapping produces large scale chaos (Chirikov, 
1960). The overlapping of resonances produces a chaotic 
band whose extent in perihelion distance away from the 
planet is proportional to the planet mass at the 2/7 power 
(Wisdom, 1980). To date, the most extended analytic cal- 
culation of the width of the mean motion resonances with 
Neptune up to order 50 has been done by D. Nesvorny, 
and the result -in good agreement with the stability bound- 
ary observed in Fig. [l]- is published electronically at 
[http://www.boulder.swri.edu/~davidn/kbminr/kbmmr.htmll 
The semi-major axis of the objects that are above the res- 
onance overlapping limit evolves by 'jumps', passing from 
the vicinity of one resonance to another, mostly during a 
conjunction with the planet. Given that the eccentricity of 
Neptune's orbit is small, the Tisserand parameter 

T=^ + 2./— (l-e2)cosz 



(aN denoting the semi-major axis of Neptune, a, e, i the 
semi-major axis, eccentricity and incUnation of the object) 
is approximately conserved. Thus, the eccentricity of the 
object's orbit has 'jumps' correlated to those of the semi- 
major axis, and the perihelion distance remains roughly 
constant. Consequently, the object wanders over the (a, e) 
plane and is effectively a member of the Scattered disk. In 
conclusion, the boundary between the black and the colored 
region in Fig. [T] marks the boundary of the Scattered disk 
and has a complicated, fractal structure, which justifies the 
use of numerical simulations in order to classify the objects 
(see Chapters by Gladman et al. and Gomes et al). 

Not all bodies with q < 35 AU are unstable, though. 
The exception is those objects deep inside low-order mean- 
motion resonances with Neptune. These objects, despite 
approaching (or even intersecting) the orbit of Neptune at 
perihelion, never get close to the planet. The reason for this 
can be understood with a little algebra. If a body is in a 
fcjv:fc resonance with Neptune, the ratio of its orbital period 
P to Neptune's Pn is, by definition, equal to k/k^. In 
this case, and assuming that the planet is on a quasi-circular 
orbit and the motions of the particle and of the planet are 
co-planar, the angle 

a = kpfX^ — kX + (k — k]\i)w, (1) 

(where A jv and A denote the mean longitudes of Neptune 
and of the object and zu is the object's longitude of peri- 
helion) has a time derivative which is zero on average, so 
that it librates around an equilibrium value, say cTstab (see 
Fig. |2]i. The radial distance from the planet's orbit is min- 
imized when the object passes close to perihelion. Perihe- 
lion passage happens when A = tn. When this occurs, from 
([T]i we see that the angular separation between the planet 
and the object, Ajv — A, is equal to a/k^. For small am- 
plitude librations, a ~ astah', because Ustab is typically far 
from see (Fig.|2]i, we conclude that close encounters can- 
not occur (Malhotra, 1996). Conversely, if the body is not 
in resonance, a circulates (Fig. |2]l. So, eventually it has to 
pass through 0, which brings the object to be in conjuction 
with Neptune during its closest approach to the planet's or- 
bit. Thus, close encounters are possible, if the object's per- 
ihelion distance is small enough. 

For most mean motion resonances, (Tstab — 180°. How- 
ever, this is not true for the resonances of type l:fc. In these 
resonances, if the eccentricity of the body is not very small, 
there are two stable equilibria at (Tstab = 180° ± 5, with 
5 ~ 60°, while a — 180° is an unstable equilibrium (Mes- 
sage, 1958; Beauge, 1994; see Fig. |2]i. Thus, bodies with 
small amplitudes of libration, necessarily librate asymmet- 
rically in (J relative to the (0, 2ti) interval. Symmetric libra- 
tions are possible only for large amplitude librators. 

A detailed exploration of the stability region inside the 
two main mean-motion resonances of the Kuiper belt, the 
2:3 and 1:2 resonances with Neptune, has been done in 
Nesvorny and Roig, 2000, 2001). In general, they found 
that orbits with large amplitude of libration and moderate to 
large eccentricities are chaotic, and eventually escape from 
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Fig. 2. — The dynamics in the 1:2 resonance with Nep- 
tune, in e cos (7, e sin (T coordinates. The motions follow 
the closed curves plotted in the figure. Librations occur on 
those curves that do not enclose the origin of the figure. 
Notice two unstable equilibrium points on the e sin cr = 
line. Each unstable equilibrium is the origin of a critical 
curve called separatrix, plotted in bold in the figure. The 
one with origin at the unstable point at cr = separates 
resonant from not resonant orbits, for which a respectively 
librates or circulates. The separatrix with origin at the un- 
stable point in (7 = 180° delimits two islands of libration 
around each of the asymmetric stable equiUbria. 

the resonance, joining the scattered disk population. Con- 
versely, orbits with small eccentricity or small libration am- 
plitude are stable over the age of the Solar System. At large 
eccentricity, only asymmetric librations are stable in the 1 :2 
resonance. 

Mean motion resonances are not the only important 
agent structuring the dynamics in the Kuiper belt. In fig. [U 
one can see that there is a dark region that extends down to 
e = when 40 < a < 42 AU. The instability in this case 
is due to the presence of the secular resonance that occurs 
when w ^ wn, where is the perihelion longitude of 
Neptune. More generally, secular resonances occur when 
the precession rate of the perihelion or of the longitude of 
the node of an object is equal to the mean precession rate of 
the perihelion or the node of one of the planets. The secular 
resonances involving the perihelion precession rates excite 
the eccentricities, while those involving the node precession 
rates excite the inclinations {Williams and Faulkner, 1981; 
Morbidelli and Henrard, 1991). 

The location of secular resonances in the Kuiper belt 
has been computed in Knezevic et al. (1991). This work 
showed that both the secular resonance with the perihe- 
lion and that with the node of Neptune are present in the 
40 < a < 42 AU region, for i < 15°. Consequently, a low- 
inclination object in this region undergoes large variations 
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Fig. 3. — The evolution of objects initially at e 0.015 
and semi-major axes distributed in the 36.5-39.5 AU range. 
The dots represent the proper semi-major axis and the ec- 
centricity of the objects -computed by averaging their a 
and e over 10 My time intervals- over the age of the So- 
lar System. They are plotted in gray after the perihe- 
lion has decreased below 32 AU for the first time. La- 
bels Nk^'-k denote the fcjv:fc two-body resonances with 
Neptune. Labels fc7vN+/c[/U+fc denote the three-body res- 
onances with Uranus and Neptune, corresponding to the 
equality k^X-^ + kuX^j + kX = 0. Reprinted from Afeivorny 
andRoig, (2001). 

in orbital eccentricity so that - even if the initial eccentricity 
is zero - the perihelion distance eventually decreases below 
35 AU, and the object enters the Scattered disk {Holman 
and Wisdom, 1993; Morbidelli et al, 1995). Conversely, a 
large inclination object in the same semi-major axis region 
is stable. Indeed, fig.[T]shows that many objects with i > 4° 
(small dots) are present between 40 and 42 AU. Only large 
dots, representing low-inclination objects, are absent. 

Another important characteristic revealed by Fig. [T| is 
the presence of narrow regions, represented by brown col- 
ored bands, where orbits become Neptune-crossing only 
after billions of years of evolution. The nature of these 
weakly unstable orbits remained mysterious for several 
years. Eventually, it was found (Nesvorny andRoig, 2001) 
that they are, in general, associated either with high or- 
der mean-motion resonances with Neptune (i.e. resonances 
for which the equivalence kX = k^XN holds only for 
large values of the integer coefficients fc, fc^v) or three-body 
resonances with Uranus and Neptune (which occur when 
kX + kpfXN + kijXjj = occurs for some integers k,kN 
and ku). 

The dynamics of objects in these resonances is chaotic 
due to the non-zero eccentricity of the planetary orbits. The 
semi-major axis of the objects remain locked at the corre- 
sponding resonant value, while the eccentricity of their or- 
bits slowly evolves. In an (a, e)-diagram like fig. |3] each 
object's evolution leaves a vertical trace. This phenomenon 
is called chaotic diffusion. Eventually the growth of the 
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eccentricity can bring the diffusing object to decrease the 
perihelion distance below 35 AU. These resonances are too 
weak to offer an effective protection against close encoun- 
ters with Neptune ((Tstah /kN is a small quantity because 
fc^v is large), unlike the low order resonances considered 
above. Thus, once the perihelion distance becomes too low, 
the encounters with Neptune start to change the semi-major 
axis of the objects, which leave their original resonance and 
evolve - from that moment on - in the Scattered disk. 

Notice from fig.[3]that some resonances are so weak that, 
despite forcing the resonant objects to diffuse chaotically, 
they cannot reach the g = 35 AU curve within the age of 
the Solar System. Therefore, although these objects are not 
stable from the dynamics point of view, they can be consider 
that way from the astronomical perspective. 

Notice also that chaotic diffusion is effective only for 
selected resonances. The vast majority of the simulated ob- 
jects are not affected by any macroscopic diffusion. They 
preserve their initial small eccentricity for the entire age of 
the Solar System. Thus, the current moderate/large eccen- 
tricities and inclinations of most of the Kuiper belt objects 
cannot be obtained from primordial circular and coplanar 
orbits by dynamical evolution in the framework of the cur- 
rent orbital configuration of the planetary system. Like- 
wise, the region beyond the 1:2 mean-motion resonance 
with Neptune is totally stable up to an eccentricity of ^ 0.3 
(Fig. [T]). As a result, the absence of bodies beyond 48 
AU cannot be explained by current dynamical instabilities. 
Therefore, these (and other) intriguing properties of the 
Kuiper belt's structure must, instead, be explained within 
the framework of the formation and primordial evolution of 
the Solar System. These topics will be treated in the next 
sections. 

3. Kuiper belt properties acquired during a primor- 
dial age 

From the current dynamical structure of the Kuiper belt, 
we conclude that the properties that require an explanation 
in the framework of the primordial Solar System evolution 
are: 

/) The existence of conspicuous populations of objects in 
the main mean motion resonances with Neptune (2:3, 3:5, 
4:7, 1:2, 2:5, etc.). The dynamical analysis presented above 
shows that these resonances are stable, but does not explain 
how and why objects populated these resonances on orbits 
with eccentricities as large as allowed by stability consider- 
ations. 

//) The excitation of the eccentricities in the classical belt, 
which we define here as the collection of non-resonant ob- 
jects with 42 < a < 48 AU and q > 37 AU. The median ec- 
centricity of the classical belt is ^ 0.07. It should be noted, 
however, that the upper eccentricity boundary of this popu- 
lation is set by the long-term orbital stability of the Kuiper 
belt (see Fig.[T]), and thus this semi-major axis region could 
have contained at some time in the past objects with much 



larger eccentricities. In any case, even if the current me- 
dian eccentricity is small, it is nevertheless much larger (an 
order of magnitude or more) than the one that must have 
existed when the KBOs formed. The current dynamics are 
stable, so that, without additional stirring mechanisms, the 
primordial small eccentricities should have been preserved 
to modern times. 

///) The peculiar (a, e) distribution of the objects in the clas- 
sical belt (see Fig. [Hi. In particular, the population of objects 
on nearly-circular orbits (e < 0.05) effectively ends at about 
44 AU, and beyond this location the eccentricity tends to 
increase with semi-major axis. If this were simply the con- 
sequence of an observational bias that favors the discovery 
of objects on orbits with smaller perihelion distances, we 
would expect that the lower bound of the a-e distribution in 
the 44^8 AU would follow a curve of constant q. This is 
not the case. Indeed, the eccentricity of this boundary grows 
more steeply with semi-major axis than this explaination 
would predict. Thus, the apparent relative under-density of 
objects at low eccentricity in the region 44 < a < 48 AU 
is likely to be a real feature of the Kuiper belt distribution. 
This under-density cannot be explained by a lack of stability 
in this region. 

iv) The outer edge of the classical belt (Fig. [T|). This edge 
appears to be precisely at the location of the 1:2 mean- 
motion resonance (MMR) with Neptune. Only large ec- 
centricity objects, typical of the scattered disk or of the de- 
tached population (see Chapter by Gladman et al. for a 
definition of these populations) seem to exist beyond this 
boundary (Fig.[T]). Again, the under density (or absence) of 
low eccentricity objects beyond the 1 :2 MMR cannot be ex- 
plained by observational biases (Trujillo and Brown, 2001; 
Allen et al., 2001, 2002; see also Chapter by Kavelaars et 
al). As Fig. [U shows, the region beyond the 1:2 resonance 
looks stable, even at moderate eccentricity. So, a primordial 
distant population should have remained there. 

v) The inclination distribution in the classical belt. The ob- 
servations (see FiglU show a clump of objects with i < 4°. 
However, there are also several objects with much larger 
inclinations, up to z ^ 30°, despite the fact that an object's 
inclination does not change much in the current Solar Sys- 
tem. Observational biases definitely enhance the low incli- 
nation clump relative to the large inclination population (the 
probability of discovery of an object in an ecliptic survey 
is roughly proportional to 1/ sin(i)). However, the clump 
persists even when the biases are taken into account. Brown 
(2001) argued that the debiased inclination distribution is 
bimodal and can be fitted with two Gaussian functions, one 
with a standard deviation cr ^ 2° for the low-inclination 
core, and the other with <j ^ 12° for the high inclination 
population (see also Chapter by Kavelaars et al). Since 
the work of Brown, the classical population with i < 4° is 
called the 'cold population', and the remaining one is called 
the 'hot population' (see Chapter by Gladman et al for 
nomenclature issues). 
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Fig. 4. — The semi-major axis — inclination distribution of 
all well observed Kuiper belt objects. The important mean 
motion resonances are also shown. 



vi) The correlations between physical properties and orbital 
distribution. The cluster of low inclination objects visible 
in the (a, i) distribution disappears if one selects only ob- 
jects with absolute magnitude H < 6 (Levison and Stern, 
2001 HI This implies that intrinsically bright objects are un- 
der represented in the cold population. Grundy et al. (2005) 
have shown that the objects of the cold population have a 
larger albedo, on average, than those of the hot population. 
Thus, the correlation found by Levison and Stern implies 
that the hot population contains bigger objects. Bernstein 
et al. (2004) showed that the hot population has a shallower 
H distribution than the cold population, which is consistent 
with the absence of the largest objects in the cold belt. In ad- 
dition, there is a well known correlation between color and 
inclination (see Chapterhy Doressoundiram et al.). The hot 
population objects show a wide range of colors, from red to 
gray. Conversely, the cold population objects are mostly 
red. In other words, the cold population shows a significant 
deficit of gray bodies relative to the hot population. The dif- 
ferences in physical properties argue that the cold and the 
hot populations have different origins. 

vii) The mass deficit of the Kuiper belt. The current mass 
of the Kuiper belt is very small. Estimates range from 
0.01 Earth masses (M®) {Bernstein et al, 2004) to 0.1 
(Gladman et al, 2001). The uncertainty is due mainly to the 
conversion from absolute magnitudes to sizes, assumptions 
about bulk density, and ambiguities in the size distribution 
(see Chapter by Petit et al.). Whatever the exact real to- 
tal mass, there appears to be a significant mass deficit (of 
2-3 orders of magnitude) with respect to what models say 



^The absolute magnitude is the brightness that the object would have if it 
were viewed at 1 AU, with the observer at the Sun. It relates to size by the 
formula LogD^ = 6.244 — O.AH — Log{p), where D is the diameter in 
kilometers and p is the albedo. 



is needed in order for the KBOs to accrete in situ. In par- 
ticular, in order to grow the objects that we see within a 
reasonable time (10^-10^ My), the Kuiper belt must have 
consisted of about 10 to 30 of solid material in a dy- 
namically cold disk (Stern, 1996; Stern and Colwell, 1997a, 
1997b; Kenyan and Luu, 1998, 1999a, 1999b; Kenyan and 
Bramley, 2004a). If most of the Kuiper belt is currently 
stable, and therefore objects do not escape, what depleted 
> 99.9% of the Kuiper belt primordial mass? 

All these issues provide us with a large number of clues 
to understand what happened in the outer Solar System dur- 
ing the primordial era. Potentially, the Kuiper belt might 
teach us more about the formation of the giant planets than 
the planets themselves. This is what makes the Kuiper belt 
so important and fascinating for planetary science. 

4. Origin of the outer edge of the Kuiper belt 

The existence of an outer edge of the Kuiper belt is very 
intriguing. Several mechanisms for its origin have been pro- 
posed, none of which have resulted yet in a general consen- 
sus among the experts in the field. These mechanisms can 
be grouped in three classes. 

Class I: Destroying the distant planetesimal disk It has 

been argued in Brunini and Melita (2002) that a Martian 
mass body residing for 1 Gy on an orbit with a ^ &{) AU 
and e ^ 0.15-0.2 could have scattered most of the Kuiper 
belt bodies originally in the 50-70 AU range into Neptune- 
crossing orbits, leaving this region strongly depleted and 
dynamically excited. It might be possible (see Chapter by 
Kavelaars et al.) that the apparent edge at 50 AU is simply 
the inner edge of such a gap in the distribution of Kuiper 
belt bodies. A main problem with this scenario is that there 
are no evident dynamical mechanisms that would ensure the 
later removal of the massive body from the system. In other 
words, the massive body should still be present, somewhere 
in the ~ 50 — 70 AU region. A Mars-size body with 4% 
albedo at 70 AU would have apparent magnitude brighter 
than 20. In addition, its inclination should be small, both in 
the scenario where it was originally a Scattered disk object 
whose eccentricity (and inclination) were damped by dy- 
namical friction (as envisioned by Brunini and Melita), and 
in the one where the body reached its required heliocen- 
tric distance by migrating through the primordially massive 
Kuiper belt (see Games et al, 2004). Thus, in view of its 
brightness and small inclination, it is unlikely that the puta- 
tive Mars-size body could have escaped detection in the nu- 
merous wide field ecliptic surveys that have been performed 
up to now, and in particular in that described in Trujilla and 
Brawn (2003). 

A second possibility for destroying the Kuiper belt be- 
yond the observed edge is that the planetesimal disk was 
truncated by a close stellar encounter The eccentricities 
and inclinations of the planetesimals resulting from a stellar 
encounter depend critically on a/D, where a is the semi- 
major axis of the planetesimal and D is the closest helio- 
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centric distance of the stellar encounter (Ida et ai, 2000; 
Kobayashi and Ida, 2001). An encounter with a solar-mass 
star at ^ 200 AU would make most of the bodies beyond 
50 AU so eccentric that they intersect the orbit of Neptune, 
which would eventually produce the observed edge (Melita 
et ai, 2002). An interesting constraint on the time at which 
such an encounter occurred is set by the existence of the 
Oort cloud. It was shown in Levison et al. (2004) that the 
encounter had to occur much earlier than 10 My after 
the formation of Uranus and Neptune, otherwise most of 
the existing Oort cloud would have been ejected to inter- 
stellar space. Moreover, many of the planetesimals in the 
Scattered disk at that time would have had their perihelion 
distance lifted beyond Neptune, decoupling them from the 
planet. As a consequence, the detached population, with 
50 < a < 100 AU and 40 < g < 50 AU, would have 
had a mass comparable or larger than that of the resulting 
Oort cloud, hardly compatible with the few detections of de- 
tached objects achieved up to now. Finally, this mechanism 
predicts an correlation between inclination and semi-major 
axis that is not seen. 

A way around the above problems could be achieved if 
the encounter occurred during the first million years of So- 
lar System history (Levison et al., 2004). At this time, the 
Sun was still in its birth cluster making such an encounter 
likely. However, the Kuiper belt objects were presumably 
not yet fully formed [Stern, 1996; Kenyan and Luu, 1998), 
and thus an edge to the belt would form at the location of 
the disk where eccentricities are ~ 0.05. Interior to this lo- 
cation collisional damping is efficient and accretion can re- 
cover from the encounter, beyond this location the objects 
rapidly grind down to dust (Kenyan and Broomley, 2002). 
If this scenario is true, the stellar passage cannot be respon- 
sible exciting the Kuiper belt because the objects that we 
observe there did not form until much later. 

According to the analysis done in Levison et al. (2004), 
an edge-forming stellar encounter should not be the respon- 
sible for the origin of the peculiar orbit of Sedna (a = 
484 AU and q — 76 AU), unlike the proposed in Kenyan and 
Broomley, (2004b). In fact, such a close encounter would 
also produce a relative overabundance of bodies with peri- 
helion distance similar to that of Sedna but with semi-major 
axes in the 50-200 AU range (Morbidelli and Levison, 
2004). These bodies have never been discovered, despite 
the fact that they should be easier to find than Sedna be- 
cause of their shorter orbital period (see Chapter by Gomes 
et al). 

Class II: Forming a bound planetesimal disk from an ex- 
tended gas-dust disk In Weidenschilling (2003), it was 
shown that the outer edge of the Kuiper belt might be the 
result of two facts: i) accretion takes longer with increas- 
ing heliocentric distance and ii) small planetesimals drift 
inwards due to gas drag. According to Weidenschilling's 
models, this leads to a steepen of the radial surface density 
gradient of solids. The edge effect is augmented because. 
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Fig. 5. — Top: the time evolution of the surface density 
of solids. Bottom: the size distribution as a function of 
heliocentric distance. From Weidenschilling (2003). 

at whatever distance large bodies can form, they capture 
the ^meter-sized bodies spiraling inwards from farther out. 
The net result of the process, as shown by numerical model- 
ing in Weidenschilling (2003; see Fig.|5]l, is the production 
of an effective edge, where both the surface density of solid 
matter and the mean size of planetesimals decrease sharply 
with increasing distance. 

A somewhat similar scenario has been proposed in 
Youdin and Shu (2002). In their model, planetesimals 
formed by gravitational instability, but only in regions of 
the Solar nebula where the local solid/gas ratio was ^ 3 
times that of the Sun (Sekiya, 1983). According to the au- 
thors, this large ratio occurs because of a radial variation 
of orbital drift speeds of millimeter-sized particles induced 
by gas drag. This drift also acts to steeping the surface 
density distribution of the disk of solids. This means that 
at some point in the nebula, the solid/gas ratio falls below 
the critical value to form planetesimals, so that the resulting 
planetesimal disk would have had a natural outer edge. 

A third possibility is that planetesimals formed only 
within a limited heliocentric distance because of the ef- 
fect of turbulence. If turbulence in protoplanetary disks 
is driven by magneto-rotational instability (MRI), one can 
expect that it was particularly strong in the vicinity of the 
Sun and at large distances (where solar and stellar radiation 
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could more easily ionize the gas), while it was weaker in 
the central, optically thick region of the nebula, known as 
the 'dead zone' (Stone et al., 1998). The accretion of plan- 
etesimals should have been inhibited by strong turbulence, 
because the latter enhanced the relative velocities of the 
grains. Consequently, the planetesimals could have formed 
only in the dead zone, with well defined outer (and inner) 
edge(s). 

Class III: IVuncating the original gas disk The detailed 
observational investigation of star formation regions has re- 
vealed the existence of many proplyds (anomalously small 
protoplanetary disks). It is believed that these disks were 
originally much larger, but in their distant regions the gas 
was photo-evaporated by highly energetic radiation emit- 
ted by the massive stars of the cluster (Adams et al, 2004). 
Thus, it has been proposed that the outer edge of the Kuiper 
belt reflects the size of the original Solar System proplyd 
(Hollenbach et al, 2004). 

There is also the theoretical possibility that the disk was 
born small and did not spread out substantially under its 
own viscous evolution (Ruden and Pollack, 1991). In this 
case, no truncation mechanism is needed. Observations, 
however, don't show many small disks, other than in clus- 
ters with massive stars. 

In all the scenarios discussed above, the location of the 
edge can be adjusted by tuning the relevant parameters of 
the corresponding model. In all cases, however, Neptune 
played no direct role in the edge formation. In this context, 
it is particularly important to remark (as seen in fig.[T]l that 
the edge of the Kuiper belt appears to coincide precisely 
with the location of the 1:2 mean-motion resonance with 
Neptune. This suggests that, whatever mechanism formed 
the edge, the planet was able to adjust the final location of 
the outer boundary through gravitational interactions. We 
will return to this in section Notice that a planetesimal 
disk truncated at ^ 34 AU has been recently postulated in 
order to explain the dust distribution in the AU Mic system 
(Augereau and Beust, 2006). 

5. The role of Neptune's migration 

It was shown in Fernandez and Ip (1984) that, in the ab- 
sence of a massive gas disk, while scattering away the pri- 
mordial planetesimals from their neighboring regions, the 
giant planets had to migrate in semi-major axis as a con- 
sequence of angular momentum conservation. Given the 
configuration of the giant planets in our Solar System, this 
migration should have had a general trend (see Levison et 
al, 2006 for a review). Uranus and Neptune have diffi- 
culty ejecting planetesimals onto hyperbolic orbits. Apart 
from the few percent of planetesimals that can be perma- 
nently stored in the Oort cloud or in the Scattered disk, the 
remaining planetesimals (the large majority) are eventually 
scattered inwards, towards Saturn and Jupiter. Thus, the ice 
giants, by reaction, have to move outwards. Jupiter, on the 
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Fig. 6. — The final distribution of Kuiper belt bodies ac- 
cording to the sweeping resonances scenario (courtesy of 
R. Malhotra). This simulation was done by numerically in- 
tegrating, over a 200 My time-span, the evolution of 800 
test particles on initially quasi-circular and coplanar orbits. 
The planets are forced to migrate by a quantity Aa (equal 
to -0.2 AU for Jupiter, 0.8 AU for Saturn, 3 AU for Uranus 
and 7 AU for Neptune) and approach their current orbits 
exponentially as a{t) = a^o — Aaexp(— i/4Mi/), where 
Ooo is the current semi-major axis. Large solid dots repre- 
sent 'surviving' particles (i.e., those that have not suffered 
any planetary close encounters during the integration time); 
small dots represent the 'removed' particles at the time of 
their close encounter with a planet (e.g. bodies that entered 
the Scattered disk and whose evolution was not followed 
further). In the lowest panel, the solid line is the histogram 
of semi-major axes of the 'surviving' particles; the dotted 
line is the initial distribution. The locations of the main 
mean-motion resonances are indicated above the top panel. 

other hand, eventually ejects from the Solar System almost 
all of the planetesimals that it encounters: thus it has to 
move inwards. The fate of Saturn is more difficult to pre- 
dict, a priori. However, numerical simulations show that 
this planet also moves outwards, although only by a few 
tenths of an AU for reasonable disk's masses (Hahn and 
Malhotra, 1999; Gomes et al, 2004). 

5.1. The resonance sweeping scenario 

In Malhotra (1993, 1995) it was realized that, follow- 
ing Neptune's migration, the mean-motion resonances with 
Neptune also migrated outwards, sweeping through the pri- 
mordial Kuiper belt until they reached their present posi- 
tions. From adiabatic theory (see for instance Henrard, 
1982), some of the Kuiper belt objects over which a mean- 
motion resonance swept, were captured into resonance; 
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they subsequently followed the resonance through its mi- 
gration, with ever increasing eccentricities. In fact, it can 
be shown {Malhotra, 1995) that, for a kN : k resonance, 
the eccentricity of an object grows as 



where is the semi-major axis that the object had when it 
was captured in resonance and a is its current semi-major 
axis. This relationship neglects secular effects inside the 
mean motion resonance, that can be important if Neptune's 
eccentricity is not zero and its precession frequencies are 
comparable to those of the resonant particles (Levison and 
Morbidelli, 2003). This model can account for the exis- 
tence of the large number of Kuiper belt objects in the 2:3 
mean-motion resonance with Neptune (and also in other 
resonances), and can explain their large eccentricities (see 
fig- 01. Assuming that all objects were captured when their 
eccentricities were close to zero the above formula indicates 
that, Neptune had to have migrated ^ 7 AU in order to re- 
produce quantitatively the observed range of eccentricities 
(up to ^ 0.3) of the resonant bodies. 

In Malhotra (1995), it was also shown that the bodies 
captured in the 2:3 resonance can acquire large inclinations, 
comparable to those of Pluto and other objects. The mecha- 
nisms that excite the inclination during the capture process 
have been investigated in detail in Gomes (2000), who con- 
cluded that, although large inclinations can be achieved, the 
resulting proportion of high inclination versus low inclina- 
tion bodies, as well as their distribution in the e-i plane, 
does not reproduce the observations well. We will return to 
this issue in section ls!2l 

The mechanism of adiabatic capture into resonance re- 
quires that Neptune's migration happened very smoothly. 
If Neptune had encountered a significant number of large 
bodies, its jerky migration would have jeopardized the cap- 
ture into resonances. For instance, direct simulations of 
Neptune's migration in Hahn and Malhotra (1999) - which 
modeled the disk with Lunar to Martian-mass planetesi- 
mals - did not result in any permanent captures. Adiabatic 
captures into resonance can be seen in numerical simula- 
tions only if the disk is modeled using many more plan- 
etesimals with smaller masses (Gomes, 2003; Gomes et ai, 
2004). The constraint set by the capture process on the max- 
imum size of the planetesimals that made up the bulk of the 
mass in the disk has been recently estimated in Murray- 
Clay and Chiang (2006). They found that resonance cap- 
ture due to Neptune's migration is efficient if the bulk of 
the disk particles was smaller than ^ 100 km and the frac- 
tion of disk mass in objects with sizes > 1000 km was less 
than a few percent. This result appears too severe, because 
the results in Gomes, (2003) and Gomes et ai, (2004) show 
that resonant captures occur in disks entirely constructed of 
Pluto-mass objects, although probably with a smaller effi- 
ciency than required in Murray-Clay and Chiang's work. 

If migration really happened smoothly, Murray-Clay 
and Chiang (2006) worked out a constraint on the migra- 



tion rate. Remember that there are two islands of libration 
in the 1:2 resonance with Neptune (see Fig. |2|l. If Nep- 
tune's migration occurs in less than 10 My, they showed 
that objects captured in the 1:2 resonance should preferen- 
tially be in the trailing island (that wehre a Ubrates around 
a value CTgtab > ti")- Converesely, most of the observed ob- 
jects are in the leading island. This, at first sight, points to 
a slow Neptune's migration. As we will see below, how- 
ever, resonant objects can be captured also from the scat- 
tered disk, and Neptune's migration might have been very 
different from what was originally envisioned. Thus, it is 
unclear which kind of constraint is provided by the internal 
distribution of the 1:2 resonant objects. 

As shown in Fig. |6] if the resonance sweeping sce- 
nario can explain the existence of the resonant popula- 
tions, it cannot explain the orbital distribution in the clas- 
sical belt, between 40 and 48 AU, nor the mass deple- 
tion of the Kuiper belt. The eccentricity excitation and, in 
particular, the inclination excitation obtained in the simu- 
lation in that region are far too small compared to those 
inferred from the observations. Thus, Hahn and Malho- 
tra (2005) suggested that resonance sweeping occurred af- 
ter that some perturbation excited, and perhaps depleted, 
the planetesimal disk. Thus, the eccentricity and inclina- 
tion distribution in the classical belt would not have been 
sculpted by the sweeping process, but would be the relic of 
such previous excitation mechanism(s). A similar conclu- 
sion was reached recently by Lykawka and Mukai (2007a), 
who found that the populations of objects in distant mean 
motion resonances with Neptune (i.e. with a > 50 AU) 
and their eccentricity-inclination-libration amplitude distri- 
butions can be explained by resonance sweeping only if the 
disk in the 40-48 AU region was already pre-excited in both 
e and i. Possible mechanisms of excitation and their prob- 
lems will be briefly discussed in sect.|6] 

5.2. The origin of the hot population 

The observation that the largest objects in the hot pop- 
ulation are bigger than those in the cold population lead 
Levison and Stern, (2001) to suggest that the hot popula- 
tion formed closer to the Sun and was transported outward 
during the final stages of planet formation. Gomes (2003) 
showed that the simple migration of Neptune described in 
the previous section could accomplish this feat. 

In particular. Gomes (2003) studied the migration of the 
giant planets through a disk represented by 10,000 parti- 
cles, a much larger number than had previously been at- 
tempted. In Gomes' simulations, during its migration Nep- 
tune scattered the planetesimals and formed a massive Scat- 
tered disk. Some of the scattered bodies decoupled from the 
planet, decreasing their eccentricities through interactions 
with some secular or mean-motion resonance (see Chapter 
by Gomes et al. for a detailed discussion of how resonances 
can decrease the eccentricities). If Neptune were not mi- 
grating, the decoupled phases would have been transient. 
In fact, the dynamics are reversible, so that the eccentricity 
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Fig. 7. — The orbital distribution in the classical belt ac- 
cording to the simulations in Gomes (2003). The dots de- 
note the population that formed locally, which is only mod- 
erately dynamically excited. The crosses denote the bodies 
that were originally inside 30 AU. Therefore, the resulting 
Kuiper belt population is the superposition of a dynamically 
cold population and a dynamically hot population, which 
gives a bimodal inclination distribution. The dotted curves 
in the eccentricity vs. semi-major axis plot correspond to 
g = 30 AU and q ^ 35 AU. 



would have eventually increased back to Neptune-crossing 
values. However, Neptune's migration broke the reversibil- 
ity, and some of the decoupled bodies managed to escape 
from the resonances and remained permanently trapped in 
the Kuiper belt. As shown in fig. |7] the current Kuiper belt 
would therefore be the result of the superposition in (a, e)- 
space of these bodies with the local population, originally 
formed beyond 30 AU. 

The local population stayed dynamically cold, in partic- 
ular in inclination, because its objects were only moderately 
excited by the resonance sweeping mechanism, as in fig.H] 
Conversely, the population captured from the scattered disk 
had a much more extended inclination distribution, for two 
reasons:f/j the inclinations got excited during the scattered 
disk phase before capture and (ii) there was a dynamical 
bias in favor of the capture of high-inclination bodies, be- 
cause at large-z the ability of mean motion resonances to de- 
crease the eccentricity is enhanced (see Chapter by Gomes 
et al). Thus, in Gomes model the current cold and hot pop- 
ulations should be identified respectively with the local pop- 
ulation and with the population trapped from the scattered 
disk. 

This scenario is appealing because, assuming that the 
bodies' color varied in the primordial disk with heliocentric 
distance, it qualitatively explains why the scattered objects 



and hot classical belt objects - which mostly come from 
regions inside ^--^30 AU - appear to have similar color dis- 
tributions, while the cold classical objects - the only ones 
that actually formed in the trans-Neptunian region - have a 
different distribution. Similarly, assuming that at the time of 
Neptune's migration the maximum size of the objects was a 
decreasing function of their initial heliocentric distance, the 
scenario also explains why the biggest Kuiper belt objects 
are all in the hot population. 

As Fig. |7] shows, there may be some quantitative prob- 
lems in the reproduction of the orbital distribution of the hot 
Kuiper belt in Gomes scenario (for instance the perihelion 
distances appear to be somewhat too low). However, an is- 
sue of principle concerns the relative weight between the 
hot and cold populations. In Gomes's simulations, only a 
fraction of a percent of the original scattered disk remained 
trapped in the hot population. On the other hand, the cold 
population was not depleted by the resonance sweeping, so 
that it retained most of the original objects. Thus, if the lo- 
cal population was similar in size distribution and number 
density to the planetesimal disk from which the scattered 
disk was extracted, it should outnumber the hot population 
by a huge factor (of ^ 1000). So, in order to obtain a final 
inclination distribution that quantitatively reproduces the 
debiased inclination distribution of Brown (2001), Gomes 
had to scale down the number of objects in the cold popu- 
lation by an appropriate factor, assuming that some mecha- 
nism, not included in the simulation, caused a decimation, 
and hence a mass depletion, of the local population. These 
mechanisms are reviewed in the next section. 

Before concluding this section, we note that the work 
by Gomes (2003) also have important implications for the 
origin of the detached population. This issue is addressed 
in detail in the Chapter by Gomes et al. and therefore we 
do not discuss it here. 

6. The mass deficit problem 

As we described in section [3] the Kuiper belt only con- 
tains roughly 0.1% of the mass that is required to grow the 
objects that we see. So, the natural question is what hap- 
pened to all that mass. We refer to this issue as the 'mass 
deficit problem.' We now review ideas that are currently in 
the literature. 

6.1. Mass removal 

Two general scenarios have been proposed for the mass 
depletion: (/) a strong dynamical excitation of the Kuiper 
belt, which caused the ejection of most of the bodies from 
the Kuiper belt to the Neptune-crossing region, and ( //) the 
collisional comminution of most of the mass of the Kuiper 
belt into dust. We start our discussion with (i). 

Because dynamics are size-independent, a dynamical de- 
pletion scenario requires that the primordial population in 
the Kuiper belt had a size distribution similar to the one that 
(currently) exists, but with a number of objects at each size 
multiplied by the ratio between the primordial mass and the 
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current mass. Remember that, in the current Solar System 
configuration, most of the Kuiper belt is stable, so that dy- 
namical erosion cannot signficantly reduce the total number 
of objects. The idea is, therefore, that some perturbation, 
that acted in the past and is no-longer at work, strongly ex- 
cited the orbital distribution of the Kuiper belt population. 
Most of the original objects acquired Neptune-crossing ec- 
centricities, so that they were subsequently eliminated by 
the scattering action of the planets. Only a small fraction 
of the original population, corresponding to the surviving 
mass fraction, remained in the Kuiper belt on excited or- 
bits like those of the observed objects. Thus, this scenario 
aims at explaining at the same time the mass depletion of 
the Kuiper belt and its orbital excitation. 

A first dynamical depletion mechanism was proposed in 
Morbidelli and Valsecchi (1997) and later revisited in Petit 
et al. (2001). This mechanisms invokes the existence of a 
planetary embryo, with mass comparable to that of Mars or 
the Earth, in the scattered disk for ^ 10^ y. Another mecha- 
nism was proposed by Nagasawa et al. (2000) and invokes 
the sweeping of secular resonances through the Kuiper belt 
during the dispersion of the primordial gas disk. 

The problem with the dynamical depletion scenario, 
which was not immediately recognized, is that the ejection 
of a massive population of objects from the Kuiper belt to 
the Neptune-crossing region would cause Neptune to mi- 
grate into the Kuiper belt. After all, this scenario invokes a 
^ 15 Af© object to remove 3> 15 of disk material and 
angular momentum must be conserved. For instance, revis- 
iting the Petit et al. (2001) work with simulations that ac- 
count for the effect of the planetesimals on the dynamics of 
the massive bodies. Gomes et al. (2004) showed that even a 
disk containing ^ 4 Af® of material between 40 and 50 AU 
drives Neptune beyond 30 AU. This is much less than the 
mass required (10-30 M^) by models of the accretion of 
Kuiper belt bodies {Stern and Colwell, 1997a; Kenyan and 
Luu, 1999b). 

The sole possibility for a viable dynamical model of 
Kuiper belt depletion is if the objects were kicked directly 
to hyperbolic or Jupiter-crossing orbits and were eliminated 
without interacting with Neptune. Only the passage of a 
star through the Kuiper belt seems to be capable of such an 
extreme excitation {Kobayashi et al., 2005). However, the 
cold Kuiper belt would not survive in this case. 

We note in passing that, even if we ignore the problem of 
Neptune's migration, massive embryos or secular resonance 
sweeping are probably not able to reproduce the inclination 
distribution observed in the Kuiper belt. Thus, these mech- 
anisms are unlikely to be an alternative to the scenario pro- 
posed in Gomes (2003) for producing the hot classical belt. 
Consequently, the idea of Hahn and Malhotra (2005) and 
Lykawka andMukai (2007a) that the classical belt acquired 
its current excitation before Neptune's migration is not sup- 
ported, so far, by an appropriate excitation mechanism. 

The collisional grinding scenario was proposed in Stern 
and Colwell (1997b) md Davis andFarinella (1997, 1998), 
and then pursued in Kenyon and Luu (1999a) and Kenyon 



and Bromley (2002, 2004a). It is reviewed in detail in the 
Chapter by Kenyon et al.. In essence, a massive Kuiper belt 
with large eccentricities and inclinations would experience 
very intense collisional grinding. Consequently, most of the 
mass originally in bodies smaller than several tens of kilo- 
meters could be comminuted into dust, and then evacuated 
by radiation pressure and Poynting-Robertson drag. This 
would lead to a substantial depletion in mass. 

To work, the collisional erosion scenario requires that 
two essential conditions be fulfilled. First, it requires a 
pecuUar primordial size distribution, such that all of the 
missing mass was contained in small, easy-to-break objects, 
while the number of large objects was essentially identical 
to that in the current population. Some models support the 
existence of such a size distribution at the end of the accre- 
tion phase {Kenyon and Luu, 1998, 1999b). However, there 
are several arguments in favor of a completely different 
size distribution in the planetesimal disk. The collisional 
formation of the Pluto-Charon binary {Canup, 2005) and 
of the 2003 ELgi family {Barkume et al., 2006), the cap- 
ture of Triton onto a satellite orbit around Neptune {Agnor 
and Hamilton, 2006), and the fact that the Eris, the largest 
known Kuiper belt object, is in the detached population 
{Brown et al., 2005), suggest that the number of big bodies 
was much larger in the past, with as many as 1 ,000 Pluto- 
sized objects {Stern, 1991). Moreover, we have seen above 
that the mechanism of Gomes (2003) for the origin of the 
hot population also requires a disk's size distribution with 
~ 1 , 000 times more large objects than currently present 
in the Kuiper belt. Finally, Charnoz and Morbidelli (2007) 
showed that, if the size distribution required for collisional 
grinding in the Kuiper belt is assumed for the entire plan- 
etesimal disk (5-50 AU), the Oort cloud and the scattered 
disk would not contain enough comet-size objects to sup- 
ply the observed fluxes of long-period and Jupiter-family 
comets: the cometesimals would have been destroyed be- 
fore being stored in the comet reservoirs (also see Stern and 
Weissman 2001). So, to fulfill all these constraints and still 
have an effective collisional grinding in the Kuiper belt, one 
has to assume that the size distribution were totally different 
in the region of the proto-planet disk swept by Neptune and 
in the region of the disk that became the Kuiper belt. This, 
a priori, seems unlikely, given the proximity between the 
two regions; however, we will come back to this in sect. [8] 

The second essential condition for substantial collisional 
grinding is that the energy of collisions is larger than the en- 
ergy required for disruption of the targets. Thus, either the 
KBOs are extremely weak (the successful simulations in 
Kenyon and Bromley (2004) had to assume a specific energy 
for disruption that is at least an order of magnitude lower 
than predicted by the hydrodynamical, in particular smooth- 
particle hydrodynamical (SPH), simulations of fragmenta- 
tion by Benz andAsphaug, 1999), or the massive primordial 
Kuiper belt had a large dynamical excitation, with e ^ 0.25 
and/or i ^ 7° (as assumed in Stern and Colwell, 1997b). 
However, if, as we argued above, the hot population was 
implanted in the Kuiper belt via the low efficiency process 
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of Gomes, (2003), then it was never very massive and would 
not have had much effect on the collisional evolution of 
the cold population. Thus, the cold population must have 
ground itself down. This is unlikely because the excitation 
of the cold-population is significantly smaller than the re- 
quired values reported above. There is the possibility that 
the collisional erosion of the cold belt was due to the high- 
velocity bombardment by projectiles in the scattered disk. 
The scattered disk was initially massive, but its dynamical 
decay was probably too fast (~ 100 My, see Duncan and 
Levison, 1997). The collisional action of the scattered disk 
onto the cold belt was included in Charnoz and MorbidelU 
(2007) but turned out to be a minor contribution. 

6.2. Pushing out the Kuiper beh 

Given the problems explained just above, an alternative 
way of solve the mass deficit problem was proposed in Lev- 
ison and MorbidelU (2003). In this scenario, the primordial 
edge of the massive protoplanetary disk was somewhere 
around 30-35 AU and the entire Kuiper belt population - 
not only the hot component as in Gomes (2003) - formed 
within this limit and was transported to its current location 
during Neptune's migration. The transport process for the 
cold population had to be different from the one found in 
Gomes (2003) for the hot population (but still work in par- 
allel with it), because the inclinations of the hot population 
were excited, while those of the cold population were not. 

In the framework of the classical migration scenario 
(Malhotra, 1995; Gomes et ai, 2004), the mechanism pro- 
posed in Levison and MorbidelU (2003) was the follow- 
ing: the cold population bodies were initially trapped in 
the 1 :2 resonance with Neptune; then, as they were trans- 
ported outwards by the resonance, they were progressively 
released due to the non-smoothness of the planetary migra- 
tion. In the standard adiabatic migration scenario (Malho- 
tra, 1995), there would be a resulting correlation between 
the eccentricity and the semi-major axis of the released bod- 
ies. However, this correlation was broken by a secular res- 
onance embedded in the 1:2 mean-motion resonance. This 
secular resonance was generated by the objects in the res- 
onance, themselves. In particular, unlike previous studies 
of migration, Levison and Morbidelli included the mass of 
the objects in the resonance, which modified the precession 
rate of Neptune's orbit. 

Simulations of this process matched the observed (a, e) 
distribution of the cold population fairly well, while the 
initially small inclinations were only very moderately per- 
turbed. In this scenario, the small mass of the current cold 
population is simply due to the fact that only a small frac- 
tion of the massive disk population was initially trapped in 
the 1:2 resonance and then released on stable non-resonant 
orbits. The final position of Neptune would simply re- 
flect the primitive truncation of the protoplanetary disk (see 
Gomes et ai, 2004 for a more detailed discussion). Most 
importantly, this model explains why the current edge of 
the Kuiper belt is at the 1:2 mean-motion resonance with 



Neptune, despite the fact that none of the mechanisms pro- 
posed for the truncation of the planetesimal disk involves 
Neptune in a direct way (see section IH. The location of 
the edge was modified by the migration of Neptune, via the 
migration of the 1 :2 resonance. 

On the flip side, the model in Levison and Morbidelli 
2003) re-opened the problem of the origin of the differ- 
ent physical properties of the cold and hot populations, be- 
cause both would have originated within 35 AU, although 
in somewhat different parts of the disk. Moreover, Lykawka 
and Mukai (2007a) showed that this model cannot repro- 
duce the low-to-moderate inclination objects in the distant 
(i.e. beyond 50 AU) high-order mean motion resonances 
with Neptune. 

7. Effects of a Dynamical Instability in the Orbits of 
Uranus and Neptune 

The models reviewed in the previous sections assume 
that Neptune migrated outward on a nearly-circular orbit. 
However, a substantially different model of the evolution of 
the giant planets has been recently proposed. This model 
-often called the 'Nice model'- reproduces, for the first 
time, the orbital architecture of the giant planet system (or- 
bital separations, eccentricities, inclinations; Tsiganis etai, 
2005) and the capture of the Trojan populations of Jupiter 
(Morbidelli et ai, 2005) and Neptune (Tsiganis et ai, 2005; 
Sheppard and Trujillo, 2006). It also naturally supplies a 
trigger for the Late Heavy Bombardment (LHB) of the ter- 
restrial planets (Gomes et ai, 2005), and quantitatively re- 
produces most of the LHB's characteristics. 

In the Nice model, the giant planets are assumed to be 
initially on nearly-circular and coplanar orbits, with orbital 
separations significantly smaller than the ones currently ob- 
served. More precisely, the giant planet system is assumed 
to lie in the region from ^ 5.5 AU to ^ 14 AU, and Sat- 
urn is assumed to be closer to Jupiter than their mutual 1:2 
MMR. A planetesimal disk is assumed to exist beyond the 
orbits of the giant planets, on orbits whose dynamical life- 
time is at least 3 My (the supposed lifetime of the gas-disk). 
The outer edge of the planetesimal disk is assumed to lie at 
~ 34 AU and the total mass is ^35A/0 (see Fig. [8^). 

With the above configuration, the planetesimals at the in- 
ner edge of the disk evolve onto Neptune-scattering orbits 
on a timescale of a few million years. Consequently, the mi- 
gration of the giant planets proceeds at very slow rate, gov- 
erned by the slow planetesimal escape rate from the disk. 
Because the planetary system would be stable in absence 
of interactions with the planetesimals, this slow migration 
continues for a long time, slightly damping out as the unsta- 
ble disk particles are removed from the system (Fig. [8]). Af- 
ter a long time, ranging from 350 My to 1 . 1 Gy in the sim- 
ulations of Gomes et al. (2005) — which is consistent with 
the timing of the LHB, approximately 650 My after planet 
formation — Jupiter and Saturn eventually cross their mu- 
tual 1:2 mean-motion resonance (Fig.[8]3). This resonance 
crossing excites their eccentricities to values slightly larger 
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than those currently observed. The small jump in Jupiter's 
and Saturn's eccentricities drives up the eccentricities of 
Uranus and Neptune, however The ice giant's orbits be- 
come chaotic and start to approach each other. Thus, a 
short phase of encounters follows the resonance-crossing 
event. Consequently, both ice giants are scattered outward, 
onto large eccentricity orbits (e ~ 0.3-0.4) that penetrate 
deeply into the disk (Fig. [8};). This destabilizes the full 
planetesimal disk and disk particles are scattered all over 
the Solar System. The eccentricities of Uranus and Neptune 
and -to a lesser extent- of Jupiter and Saturn, are damped 
on a timescale of a few My due to the dynamical friction 
exerted by the planetesimals. Thus, the planets decouple 
from each other, and the phase of mutual encounters rapidly 
ends. During and after the eccentricity damping phase, the 
giant planets continue their radial migration, and eventually 
reach final orbits when most of the disk has been eUminated 
(Fig.®). 

The temporary large eccentricity phase of Neptune 
opens a new degree of freedom for explaining the orbital 
structure of the Kuiper belt. The new key feature to the 
dynamics is that, when Neptune's orbit is eccentric, the full 
(a, e) region up to the location of the 1 :2 resonance with the 
planet is chaotic, even for small eccentricities. This allows 
us to envision the following scenario. We assume, in agree- 
ment with several of the simulations of the Nice model, that 
the large eccentricity phase of Neptune is achieved when 
the planet has a semi-major axis of ^ 28 AU, after its last 
encounter with Uranus. In this case, a large portion of the 
current Kuiper belt is already interior to the location of the 
1:2 resonance with Neptune. Thus, it is unstable, and can 
be invaded by objects coming from within the outer bound- 
ary of the disk (i.e. ^4 AU). When the eccentricity of 
Neptune damps out, the mechanism for the onset of chaos 
disappears. The Kuiper belt becomes stable, and the objects 
that happen to be there at that time remain trapped for the 
eternity. Given that the invasion of the particles is fast and 
the damping of Neptune's eccentricity is also rapid, there 
is probably not enough time to excite significantly the par- 
ticles' orbital inclinations if Neptune's inclination is also 
small. Therefore, we expect that this mechanism may be 
able to explain the observed cold population. The hot pop- 
ulation is then captured later, when Neptune is migrating 
up to its final orbit on a low-eccentricity orbit, as in Gomes 
(2003). 

The numerical simulations of this process have been pre- 
sented in Levison et al. (2007). Fig.|9]compares with the ob- 
servations the semi-major axis vs. eccentricity distribution 
resulting from one of the simulations, 1 Gy after the giant 
planet instability. The population of quasi-circular objects 
at low inclination extends to ~ 45 AU, in nice agreement 
with the observations. The deficit of low eccentricity ob- 
jects between 45 and 48 AU is reproduced, and the outer 
edge of the classical belt is at the final location of the 1:2 
MMR with Neptune as observed. Moreover, the real Kuiper 
belt shows a population of objects with q ^ 40 AU beyond 
the 1 :2 resonance with Neptune, which are known to be sta- 
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Fig. 8. — Solar System evolution in the Nice model, (a): at 
a time close to the beginning of the evolution. The orbits of 
the giant planets (concentric circles) are very close to each 
other and are quasi-circular. They are surrounded by a disk 
of planetesimals, whose inner edge is due to the perturba- 
tions from the planets and the outer edge is assumed to be at 
34 AU. (b): immediately before the great instability. Saturn 
is about crossing the 1:2 resonance with Jupiter, (c): at the 
time of the instability. Notice that the orbits of the planets 
have become eccentric and now penetrate the planetesimal 
disk, (d): after the LHB. The planets are parked on orbits 
very similar (in terms of separation, eccentricity, and incli- 
nation) to their current ones. The massive planetesimal disk 
has been destroyed. Only a small fraction of the planetesi- 
mals remain in the system on orbits typical of the scattered 
disk, Kuiper belt, and other small body reservoirs. From 
Gomes et al. (2005). 



ble (Emelyanenko et al. 2003). This population has been 
known by several names in the literature, which include the 
fossilized scattered disk, the extended scattered disk, and 
the detached population (see Chapter by Gladman et al). 
This model reproduces this population quite well. 

Three main differences are also noticeable, though. (/) 
All the mean motion resonances are overpopulated relative 
to the classical belt. This is probably the consequence of the 
fact that in the simulations the migration of Neptune's orbit 
and its eccentricity damping were forced smoothly, through 
fake analytic terms of the equations of motion. As we said 
above, a migration with some stochastic component (due 
to the encounters with massive objects in the disk) would 
have produced fewer surviving bodies in the resonances. 
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Fig. 9. — The distribution of semi major axes and eccen- 
tricities in the Kuiper belt. Left panel: result of a simulation 
based on the Nice model. Right panel: the observed distri- 
bution (3 oppositions objects only). The vertical solid lines 
mark the main resonance with Neptune. The dotted curve 
denotes perihelion distance equal to 30 AU and the dashed 
curve delimits the region above which only high inclination 
objects or resonant objects can be stable over the age of the 
Solai- System. 
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Fig. 10. — The cumulative inclination distribution of the 
observed classical belt objects (solid curve) and that ex- 
pected from the result of our simulation, once the obser- 
vational biases are taken into account (dotted curve). 



(ii) The region above the long-dashed curve is overpopu- 
lated in the simulation. The curve represents approximately 
the boundary between the stable (yellow) and the unstable 
(black) regions in Fig.[T] Thus, if the final orbits of the gi- 
ant planets were exactly the same as the real ones and the 
simulations were extended for the age of the Solar System, 
most of the population above the curve would be depleted as 
a consequence of chaotic dynamics. (Hi) The cold Kuiper 
belt has eccentricities that are slightly too large. The me- 
dian eccentricity of the real objects with 42 < a < 48 AU 
and g > 37 AU is 0.07, while the model produces a value of 
0.10. 

Figure [TO] shows the cumulative inclination distribution 
of objects trapped in the classical belt at the end of a sim- 
ulation and compares it with the observed distribution. For 
the comparison to be meaningful, the simulated distribution 
was run through a survey bias-calculator, following the ap- 
proach of Brown (2001). The two curves are very similar 
Indeed, they are almost indistinguishable for inclinations 
less then 6 degrees. This means that the cold population, 
and the inclination distribution within it, has been correctly 
reproduced, as well as the distribution in the lower part of 
the hot population (i.e. that with 4° < « < 10°). We re- 
mark, however, a dearth of large inclination objects. This 
deficit is intriguing and unexplained, in particular given that 
the raw simulations of the Nice model (namely those in 
which the planets are not forced to migrate, but are let free 
to respond to their interactions with massive planetesimals) 
produce objects captured in the classical belt or in the de- 
tached population with inclinations up to 50° (see Fig. 3 in 
the Chapter by Gomes et ai). In general, we would ex- 



pect an inclination distribution in the hot population that is 
equivalent to that of Gomes (2003), or even more excited. 
In fact, as pointed out in Lykawka and Mukai (2007b), the 
inclinations in the scattered disk, from which the hot pop- 
ulation is derived, are restricted to be less than ^ 40° by 
the conservation of the Tisserand parameter with respect to 
Neptune, which holds if the planet is on a quasi-circular 
orbit as in the simulations of Gomes (2003). In the Nice 
model, the eccentricity of Neptune breaks the conservation 
of the Tisserand parameter, and hence, in principle, inclina- 
tions can be larger. 

The results of the simulations based on the Nice model 
also provide a qualitative explanation for the observed cor- 
relations between inclination and physical properties. The 
particles that are trapped in the cold classical belt come, al- 
most exclusively, from the outermost parts of the planetesi- 
mal disk — in particular beyond 29 AU. Conversely, a sig- 
nificant fraction of those trapped in the hot population come 
from the inner disk. Thus, if one assumes that the largest 
objects could form only in the inner part of the disk,then 
these objects can only (or predominantly) be found in the 
hot population. Similarly, if one assumes that (for some un- 
known reason) the objects from the outer part of the disk are 
red and those from the inner part are gray, the cold popula- 
tion would be composed almost exclusively of red objects, 
whereas the hot population would contain a mixture of red 
and gray bodies. 

The simulations in Levison et al. (2007) show that 50 to 
130 particles out of 60,000 are trapped in the classical belt 
(cold and hot populations together in roughly equal propor- 
tion). According to the Nice model, the original planetes- 
imal disk contained 35 M^, thus this model predicts that 
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Fig. 11. — The eccentricity-inclination distribution of the 
Plutinos. Left panel: the simulated distribution; black dots 
refer to particles from the outer disk and gray triangles to 
particles from the inner disk. Right panel: the observed 
distribution. The relative deficit of observed Plutinos at low 
eccentricity with respect to the model is probably due to ob- 
servational biases and to the criterion used to select 'reso- 
nant objects' from the simulation (see Levison et al., 2007). 



the classical Kuiper belt should currently contain between 
^ 0.02 and ~ O.OSAfg, in good agreement from observa- 
tional estimates. Of course, to be viable, the model needs 
to explain not only the total mass of the belt, but also the 
total number of bright, detectable bodies. It does this quite 
nicely if one assumes that the original disk size-distribution 
is similar to the one currently observed. As we explained in 
section l6n this is consistent with other constraints like the 
formation of the Pluto-Charon binary. Thus, the Nice model 
explains, for the first time, the mass deficit of the Kuiper 
belt and the ratio between the hot and the cold population, 
in the framework of an initial planetesimal size distribution 
that fulfills all the constraints enumerated in section |6] 

Finally, the Nice model reproduces in a satisfactory way 
the orbital distributions of the populations in the main mean 
motion resonances with Neptune. Fig. [TT] compares the 
(e, i) distribution of the Plutinos obtained in one of Levison 
et al. (2007ys simulations, against the observed distribu- 
tion. The overall agreement is quite good. In particular, this 
is the first model that does not produce an overabundance 
of resonant objects with low inclinations. 

Moreover, the left panel of Fig, [TT] uses different sym- 
bols to indicate the particles captured from the inner (a < 
29 AU) and the outer (a > 29 AU) parts of the disk. As 
one sees, the particles are very well mixed, which is in 
agreement with the absence of correlations between col- 
ors and inclinations among the Plutinos. Conversely, a 
very strong correlation was expected in the original Gomes 
(2003) scenario because a large number of low-inclination 
bodies were captured from the cold disk. The reason that 
the Nice model is so much more successful than previous 
models is that the 2:3 resonance cannot capture any objects 



via the mechanism of Malhotra (1995). This is due to the 
fact that the resonance is already beyond the disk's outer 
edge at the beginning of the simulation (i.e. after the last 
encounter of Neptune with Uranus). We believe that this 
success strongly supports that idea of a planetesimal disk 
truncated at 30-35 AU and of a 'jump' of Neptune towards 
the outer edge of the disk. 

As for the higher order resonances beyond 50 AU, the 
simulations of Levison et al. (2007) produce populations 
with moderate libration amplitudes and inclinations consis- 
tent with observations, thus satisfying the constraint posed 
by Lykawka and Mukai (2007a). 

In summary, the strength of the Nice model is that it is 
able to explain all the intriguing properties of the Kuiper 
belt, at least at a semi-quantitative level, in the framework of 
a single, unique event. That the same scenario also explains 
the orbital architecture of the giant planets, the Trojans of 
both Jupiter and Neptune, and the LHB is, of course, a non- 
negligible additional plus that should give credence to the 
model. 

Of course, the Nice model is not perfect. As we have 
seen above, the simulations thus far performed have not 
been able to simultaneously produce: 1) the very highest 
inclinations that we see, and 2) a cold belt that is cold 
enough. The numerical simulations contain some simpli- 
fications which might effect the results. In particular, mu- 
tual collisions and collective gravitational effects among the 
planetesimals are neglected. Moreover, as discussed in Lev- 
ison et al (2007), only a subset of the evolutions of the gi- 
ant planets observed in the simulations of the Tsiganis et al. 
(2005) will produce a cold classical belt. Some experiments 
overly excite inclinations so that a cold belt is not formed, 
despite producing good final planetary orbits. Nonetheless, 
we feel that the Nice model's strengths outweigh it weak- 
nesses, particularly given that other models of Kuiper belt 
formation have had much more limited success at reproduc- 
ing the observations. 

7.1. Other planetary instability models 

The Nice model is not the first model to make use of a 
temporary dynamical instability of the giant planet system 
(and probably also not the last one!). 

Thommes et al. (1999) proposed that Uranus and Nep- 
tune formed in between the orbits of Jupiter and Saturn. 
They were subsequently destabilized and scattered onto or- 
bits with larger semi-major axis and eccentricities. The dy- 
namical excitation was eventually damped by the dynami- 
cal friction exerted by a massive planetesimal disk, and the 
planets achieved stable orbits. Their simulations showed 
an interesting sculpting of the (a, e) distribution in the re- 
gion corresponding to the classical Kuiper belt. However, 
the planetesimal disk was extended to 60 AU. Thus, no 
outer edge was produced at the 1:2 resonance with Nep- 
tune and there was not enough mass depletion in the Kuiper 
belt. Moreover, since Uranus and Neptune started between 
Jupiter and Saturn they suffered much stronger encounters 
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with the gas giants than occurred in the Nice model. As a re- 
sult, the planetesimal disk needed to be much more massive 
— so massive that if the simulations had been run to com- 
pletion, Neptune would have migrated well beyond 30 AU. 

Chiang et al. (2006) have recently speculated on a sce- 
nario based on recent work by Goldreich et al. (2004a,b), 
who, from analytic considerations, predicted the formation 
of 5 planets between 20 and 40 AU. These planets remained 
stable during their formation because their orbits were con- 
tinuously damped by the dynamical friction exerted by a 
disk of planetesimals that contained more mass than the 
planetary system. These planetesimals were very small 
(sub-meter in size) and thus remain dynamically cold due 
to collisional damping. When the planets reached Neptune- 
mass, the mass of the planets and the mass of the disk be- 
came comparable, so that the planets became unstable. Gol- 
dreich et al. conjectured that 3 of the 5 planets were ejected 
and the two remaining ones stabilized on orbits comparable 
to those of Uranus and Neptune. 

Chiang et al. suggested that at the time of the instabil- 
ity, the disk contained two populations: one made up of 
~ 100 km objects and one consisting of sub-meter objects. 
The current Kuiper belt structure would be the result of the 
orbital excitation suffered during the multi-planet instabil- 
ity. The hot population would be made up of the larger ob- 
jects which were permanently excited during the instability. 
In contrast the smallest planetesimals would suffer a sig- 
nificant amount of collisional damping, which would have 
led to the eventual accretion of the cold population. Nu- 
merical simulations made by Levison and Morbidelli (2007) 
with a new code that accounts for a planetesimal disk with 
strong internal collisional damping, invalidate the Goldre- 
ich et al. proposal. It is found that a system of 5 unstable 
Neptune-mass planets systematically leads to a system with 
more than 2 planets, spread in semi-major axis well beyond 
30 AU. Thus, the architecture of the solar planetary system 
is inconsistent with Goldreich et al.'s idea. 

8. Conclusions and Discussion 

In this chapter we have tried to understand which kind 
of solar system evolution could have produced the most im- 
portant properties of the Kuiper belt: its mass deficit, its 
outer edge, the co-existence of a cold and a hot classical 
population with different physical properties, and the pres- 
ence of resonant populations. We have proceeded by basic 
steps, trying to narrow the number of possibilities by con- 
sidering one Kuiper belt feature after the other, and starting 
from the most accepted dynamical process (planet migra- 
tion) and eventually ending with a more extravagant one (a 
temporary instabihty of the giant planets). 

We have converged to a basic scenario with three ingre- 
dients: the planetesimal disk was truncated close to 30 AU 
and the Kuiper belt was initially empty; the size distribu- 
tion in the planetesimal disk was similar to the current one 
in the Kuiper belt, but the number of objects at each size 
was larger by a factor of ~ 1, 000; the Kuiper belt objects 



are just a very small fraction of the original planetesimal 
disk population and were implanted onto their current orbit 
from the disk during the evolution of the planets. A tem- 
porary high-eccentricity phase of Neptune, when the planet 
was already at ^ 28 AU -as in the Nice model- seems to 
be the best way to implant the cold population. 

The mass deficit problem is the main issue that drove 
us to this conclusion. In particular, we started from the 
consideration that, if a massive disk had extended into the 
Kuiper belt, neither collisional grinding nor dynamical ejec- 
tion could have depleted its mass to current levels. Dynam- 
ical ejection seems to be excluded by the constraint that 
Neptune did not migrate past 30 AU. Collisional grinding 
seems to be excluded by the arguments that the size distri- 
bution was about the same everywhere in the disk and that 
~ 1, 000 Pluto-size bodies had to exist in the planetary re- 
gion. 

Is there a flaw in this reasoning? Are we really sure that 
the cold population did not form in situ? The argument that 
the size distribution of the in-situ population should be sim- 
ilar to that in the region spanned by Neptune's migration ne- 
glects possible effects due to the presence of an edge. After 
all, an edge is a big discontinuity in the size and mass dis- 
tribution, so that it may not be unreasonable that the region 
adjacent to the edge had very different properties from the 
region further away from the edge. 

Our view of the Kuiper belt evolution could radically 
change if a model of accretion were developed that pro- 
duces a disk of planetesimals with a size distribution that 
changes drastically with distance, such that (a) beyond 45 
AU all objects are too small to be detected by telescope 
surveys, (b) in the 35-45 AU region the distribution of the 
largest objects is similar to that observed in the current cold 
population while most of the mass is contained in small 
bodies, and (c) within 35 AU most of the mass is contained 
in large bodies and the size distribution culminating with 
~ 1, 000 Pluto-sized objects. If this were the case, the disk 
beyond 35 AU could lose most of its mass by collisional 
grinding before the beginning of Neptune's migration, par- 
ticularly if the latter was triggered late as in our LHB sce- 
nario (see sect. 3.4 of the chapter by Kenyon et al.). Within 
35 AU, because of the different size distribution, colhsional 
grinding would have been ineffective {Charnoz and Mor- 
bidelli, 2007). Therefore, at the time of the LHB, the sys- 
tem would have been similar to the one required by the Nice 
model, in that Neptune would have seen an effective edge 
in the planetesimal mass distribution that would have kept 
it from migrating beyond 30 AU. 

Whether the spatial variation we described above of the 
size distribution in the planetesimal disk is reasonable or 
not is beyond our current understanding. The chapter by 
Kenyon et al. in this book, nicely shows that the coagula- 
tion/erosion process is always on the edge of an instabil- 
ity. In fact, the dispersion velocity of the small bodies is 
of the order of the escape velocity from the largest bod- 
ies. Depending on the details of the collisional cascade 
(see sect. 3.5 in that chapter) the dispersion velocity can be 
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slightly smaller than the escape velocities (favoring the ac- 
cretion of a large number of massive bodies and producing 
a top-heavy size distribution which does not allow an effec- 
tive collisional grinding), or can be slightly larger (stalling 
runaway accretion and leaving most of the mass in small, 
easy to break bodies). Perhaps, the inner part of the disk 
was in the first regime and the outer part was in the second 
one, with a relatively sharp transition zone between the two 
parts. More work is needed to clarify the situation, with 
a close collaboration between experts of accretion and of 
dynamical evolution. 

From a purely dynamical point of view, the Nice model 
is not inconsistent with the existence of a local, low mass 
Kuiper belt population, extended up to 44^5 AU. In fact, 
Fig.[T2]compares the (a, e) distribution observed in the clas- 
sical belt with q > 38 AU (bottom panel) with the one that 
the Nice model predicts assuming that the disk was trun- 
cated at 34 AU (top panel, enlargment of the left panel in 
Fig. lU or assuming that the disk was truncated at 44 AU 
(middle panel). The two model distributions are statistically 
equivalent, and are both very similar to the observed distri- 
bution. In the case where the outer edge was placed at 44- 
45 AU about 7% of the particles initially in the Kuiper belt 
(a > 40 AU, g > 38 AU) remain there, although their or- 
bits have been modified. The others escape to larger eccen- 
tricities during the phase when the Kuiper belt is globally 
unstable, discussed in sect. 7. If 90% of the local mass es- 
capes, the local belt had to have been signficantly depleted 
before the time of the LHB, probably accounting for only 
a few tens of an Earth mass, otherwise, presumably (see 
sect. 6. 1), Neptune would have been driven past 30 AU. 

We end this chapter by encouraging observers to attempt 
to probe the regions beyond 50 AU. In particular, even if the 
absence of a population of objects similar to that in the 40- 
50 AU region is now secure, nothing is known on the pos- 
sible existence of small objects. Recent claims on a large 
number of stellar occultation events by distant ~ 100m ob- 
jects (Rogues et al., 2006; Chang et al., 2006; Georgevits, 
2006) may suggest, if confirmed (see Jones et al. 2006 
for a rebuttal of Chang et al. results), the existence of an 
extended disk of small planetesimals that did not grow to 
directly detectable sizes. If true, this would give us ex- 
tremely valuable information on the origin of the Kuiper 
belt edge. The model by Weidenschilling predicts the ex- 
istence of an extended disk of exclusively small bodies, 
although of sizes up to ^ Im rather than ^ 100 m (see 
Fig. Is). Conversely, other models on edge formation, such 
as the disk stripping by passing stars, photo-evaporation or 
turbulent stirring, would predict the total absence of objects, 
of any size. Therefore we encourage the pursuit of stellar 
occultation programs, until the real situation is clarified. 

So far, the Kuiper belt has taught a lot to us about the 
evolution of the planets. Planet migration, for instance, had 
been totally overlooked despite the pioneer work by Fer- 
nandez and Ip, until the Kuiper belt was discovered. More- 
over, as we mentioned above, we can find in the Kuiper 
belt further evidence that the giant planets passed through a 
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Fig. 12. — The semi-major axis vs. eccentricity distribution 
in the classical belt (q > 38 AU). Top panel: the distribu- 
tion of the objects that originated within 34 AU, implanted 
in the Kuiper belt at the time of the LHB. Middle panel: the 
distribution of the bodies assuming the disk originally ex- 
tended to 44 AU. Bottom panel: the observed distribution. 



temporary phase of violent instability. But, the Kuiper belt 
can do more. It can potentially teach us about planetesimal 
formation and the growth of larger objects because the ob- 
jects that inhabit it most likely probe different regions of the 
proto-planetary disk where accretion proceeded in very dif- 
ferent ways. Thus, it is a dreamed-of laboratory to test and 
calibrate the accretion models. This is probably the main 
venue for the future. 
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